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Abstract 

The dynamics of the convergence to the critical attractor for the lo- 
gistic map is investigated. At the border of chaos, when the Liapunov 
exponent is zero, the use of the non-extensive statistical mechanics for- 
malism allows to define a weak sensitivity or insensitivity to initial con- 
ditions. Using this formalism we analyse how a set of initial conditions 
spread all over the phase space converges to the critical attractor in the 
case of infinite bifurcation and tangent bifurcation points. We show that 
the phenomena is governed in both cases by a power-law regime but the 
critical exponents depend on the type of bifurcation and may also depend 
on the numerical experiment set-up. Differences and similarities between 
the two cases are also discussed. 



1 Introduction 

Nonlinear dissipative maps are widely used in the description of a variety of 
systems exhibiting complex behaviour, given their suitability for numerical sim- 
ulations. Sensitivity to initial conditions is easy to investigate in association 
to the dynamics onto the fractal attractor. Physical properties like entropies 
and relaxation phenomena are easily computed and generalised to more sophis- 
ticated systems. The logistic map is the most famous and used among the 
nonlinear maps. Nevertheless it still presents some aspect of the dynamics not 
well understood or exploited. This is the case for special points of the map, like 
the border or onset of chaos, constituted by the infinite bifurcation points or 
the tangent bifurcation points, often indicated as critical points, where the Lia- 
punov exponent vanishes and sensitivity to initial conditions is not exponential 
but polynomial in time. A vast literature [C. Tsallis ct.al., 1997; F. Baldovin 
et al., 2002; U. Tirnakli ct al., 1999; M. Coraddu et al., 2004] has recently 
tried to investigate these critical points especially in association to entropic or 
relaxation properties. The main feature of the logistic map at these points is 
a power-law sensitivity or insensitivity to initial conditions [P. Grassberger and 
M. Schcuncrt, 1981; C. Tsallis, 1988; F.A.B.F. Dc Moura et al., 2000], with con- 
sequent peculiar behaviour for the time evolution and the convergence to the 
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attractor of an ensemble of points. Recently a power-law sensitivity to initial 
conditions has been rewritten [C. Tsallis et al., 1997; M.L. Lyra and C. Tsallis 
1998] in the form 

f(t) = [i + (i-g)V] 1/(1 - ,) , (i) 

where X q defines the time scale after which the power-law behaviour begins. 
When used to investigate the expansion (weak sensitivity, namely l/(l-q)>0) 
or contraction (weak insensitivity, conversely 1 / (1-q) < 0) of a compact ensemble 
of points or grouped initial conditions this relation is connected to the evolution 
of the entropy of such a system. To be more clear let consider a partition of the 
phase space of the logistic map in Wb ox cells, where we use the parameterisation 

x n +i = l-ax 2 n (2) 

so that the phase space is— 1 < a; < +1. One of the standard ways to define 

the entropy of a system is the Boltzmann-Gibbs (BG) formulation 

w boai 

S BG = - ^Pilogfe) (3) 
i=i 

where pi defines the probability of occupation of the i-th cell as the ratio of points 
inside the cell with respect to the total number of points in the phase space. 
When the parameter a assumes values such that the map behaves chaotically, the 
Liapunov exponent is positive and there is strong sensitivity to initial conditions. 
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If one groups a set of points close together they will expand exponentially in 
time (at least for small time scales) and correspondingly the BG entropy will 
increase. It has been demonstrated [A.N. Kolmogorov, 1958; Ya. G. Sinai, 1959] 
that the rate of increase of the BG entropy, known as Kolmogorov-Sinai (KS) 
entropy and defined as 

_ .. S(t)-S(0) 

S KS = hm w (4) 

t— >oo t 

equals the Liapunov exponent under certain conditions, giving the well known 
Pcsin equality A = Sks [Ya. Pesin, 1977], so that the entropy increases linearly 
with time. Conversely, if the map behaves periodically, the Liapunov exponent 
is negative and the BG entropy will decrease. 

When the parameter a is set in a critical point the Liapunov exponent van- 
ishes and a set of different features appears. Both at the onset of chaos and at 
tangent bifurcation points new definitions of the sensitivity and of the entropies 
are suitable to describe the weak sensitive/insensitive behaviour. The term 
weak sensitivity (insensitivity) designs a separation that increases (decreases) 
not exponentially (the Liapunov exponent is zero) but polynomially in time. 

Two different numerical experiments present a notable interest with regard 
to the way the points evolve in time approaching the attractor. 

On the one extreme, it is interesting to choose (randomly or not), all the 
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initial conditions in a single cell, so that the initial entropy is zero [E. Borges et 
al., 2002; M. Coraddu et al., 2004]. The number of cells W(t) occupied at time 
t then grows as: 



where K q is the generalised KS entropy, defined as the rate of increase of the 
non-extensive generalised Tsallis entropy 



appropriate to describe the behaviour of the logistic map at critical points [C. 
Tsallis 1988]. 

On the other extreme one can choose to spread the initial conditions all over 
the whole phase space, occupying most of (or all, if the set of initial points is at 
least equal to the number of cells) the partition cells [F.A.B.F. De Moura et al., 
2000] . The fact that the parameter a has been set in a critical point allows the 
system to exhibit peculiar behaviour. The main interest of this paper consists in 
the numerical analysis of the time evolution of a set of initial conditions spread 
in the whole phase space and the way they converge to the critical attractor. 

We analyse the case of weak insensitivity and make a comparison with the 
weak sensitivity case. Due to the criticality of the system, the temporal evolu- 
tion may depend on the particular choice of the initial conditions. In our case, 



W{t) = [\ + (\-q)K q t]W-*> 



(5) 




(6) 
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with initial conditions spread over the phase space, one expects a shrink toward 
the attractor in both cases, namely a corresponding to the infinite bifurcation 
point (aoo = 1.401155189) and to the tangent bifurcation point (at g = 1.75). 

2 The Convergence to the Critical Attractor 

We analyse the temporal evolution of an ensemble of N initial conditions ran- 
domly chosen and uniformly distributed in the whole phase space. Following 
the work done in [F.A.B.F. de Moura et al.] we partitioned the phase space 
in Wbox equal cells and examined the number of occupied cells W(t) at time t 
for different ratios r = N/Wb ox . Physically r is an index giving the sampling 
ratio of the phase space. We used Wbox ranging from 5000 to 128000 cells for 
the partition and ratios r varying from r — 0.1 to r — 10000, for both cases 
o. = aoc = 1.401155189 and a — a tg — 1.75. When the convergence to the 
critical attractor is analysed using the non-extensive Tsallis entropy © , there 
is a particular value of the entropic index q [G.F.J. Ananos et al., 2004 and 
references therein] such that the entropy evolves in time at constant rate 
K q , with 

K q = km [SS) ~ S q Q0)]/t (7) 

N— >oo 

In the case of relaxation onto the attractor, the number of occupied cells de- 
creases in time giving a negative K q . In the simplest case of equiprobability, so 
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that each cell contains the same number of points, the relation 10 gives 

and the number of occupied cells is expected to present a power-law decay 

W{t) = [^(O) 1 " 9 + (1 - q)K q t] 1 / ( - 1 -^ (9) 
with asymptotic exponent a — l/(l-q) < 0. 

In order to exploit most of the features of the convergence to the critical at- 
tractor we set up a series of different numerical experiments varying W(, ox at 
constant r and vice- versa, for the two cases a = eioo and a = at g . The nu- 
merical experiments consist in selecting, using a uniform random distribution, 
at time zero, N initial points in the phase space of the logistic map x 6 [-1.1], 
that has been partitioned in Wb ox boxes. As each of the N points is iterated in 
time according to the logistic map W(t) boxes, among all the Wbox, will be not 
empty at time t. The number W(t) generally will decrease iterating the map 
and is the quantity used in the following to characterize the different regimes of 
convergence to the attractor of the N initial points. 

As a first case we analyse the infinite bifurcation point. Here the expected 
power-law behaviour (Eq. (9)) [P. Grassberger and M. Scheunert, 1981] sets up 
for a — even at low r values (r = 0.1); increasing r a superimposed log- log 
oscillation appears clearly and becomes more evident at higher rvalues (Fig. 1). 
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In these cases there is a transient behaviour, preceeding the power-law regime , 
which is longer for larger r. 

It is important to note that, in a logarithmic-to-logarithmic scale, the tran- 
sient time T crossover after which the convergence starts is translated by almost 
the same amount when r increases of a factor 10. This feature seems to persist 
up to the values of r = 10000, reached in our experiments. Higher values for 
the sampling ratio are in order to be used to verify the power-law regime in the 
limit r — > oo. 

The dependence of T crossover on the samping ratio r is clearly evidenced in 
Fig. 1. Previous works analysed values of sampling ratio up to r = 10, obtaining 
as a result a crossover time T crossover independent on r and related to the Eq. 
(??) by T crossover - wio)"- 1 ■ Ana ly sm g sampling ratios up to r = 10000, 
for many different Wbox, this picture breaks down revealing a net and regular 
variation of T crossover with r. 

In our analysis the power-law exponent seems to be not independent on r. 
The relation 

holds in the hypothesis of log-log oscillations [D. Sornette et al., 1998; C. Tsallis 
et al., 1997], where P() is a period one function and A is a scaling factor between 
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two oscillation periods. If the exponent \i is constant a plot of W(t)/t~ M should 
show a periodic function in log(t) with constant average. Figure 2 shows that the 
log-periodic oscillations exist but the function P() can be globally increasing or 
decreasing in average, depending on the value of r. This trend can be explained 
if the exponent fj, is not perfectly constant but slightly varying with r (at least 
for the considered values; it may be possible that increasing r a limit value for 
\i shows up. This value should probably be less than -0.71, around -0.8). In 
the figure we used the value \i = -0.71 according to [F.A.B.F. Dc Moura et al., 
2000] . However to the same sampling ratio r corresponds the same exponent 
irrespectively of Wb ox . 

In the case of the tangent bifurcation point the regime of convergence to 
the critical attractor also presents a power-law behaviour, but with different 
exponents (Fig. 3). Fitting the curves we found fi = -1.55 for r = 1000 and 
/j, = -1.45 for r = 10000 (Fig. 4a). This difference is not surprising, since the 
exponent is related to the structure of the chaotic attractor [F.A.B.F. De Moura 
et al., 2000]. Nevertheless a power- law behaviour was not obvious at the tangent 
bifurcation point because the critical attractor has only three fixed point and 
does not posses a multifractal structure. To our knowledge this is the first 
study of such a power-law behaviour during the convergence to the attractor 
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of a statistical ensemble of initial conditions in a critical point with not infinite 
bifurcations or multifractal structure. We found (Fig. 3), as in the case of 
the infinite bifurcation point (Fig. 1), small irregular oscillations perturbing the 
powerdaw behaviour (to not be confused with the logdog oscillations) due to the 
finite size partition of the phase space. The larger is the number of cells Wbox 
(the finer is the partition) the smaller are the oscillations around the powerdaw 
behaviour. These oscillations are irregular and grid (and/or r) depending, and 
are almost completely reduced when Wbox is big enough. 

No logdog oscillations are observed in this case (Fig. 4b), reflecting the 
different structure of the attractor. 

Again, increasing r, in the log-todog scale, the transient time T crossover after 
which the powerdaw regime of convergence begins is translated by approxima- 
tively the same amount when r is increased of a factor 10 (Fig. 3). The ratio 
W(t)/W(0) before T crossover is determined in both cases by the portion of phase 
space occupied by the attractor at the edge of chaos. 

We repeated the experiments with a different kind of initial set-up, choosing 
all the initial N points concentrated in a single cell, for a = and a — a tg , 
using as a test case a grid of 64000 cells. At the infinite bifurcation point we 
found that, after a transient time, the system reaches a powerdaw regime with a 
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different exponent with respect to the experiment with all the initial conditions 
spread in x £ [-1,1] (Fig. 5a), so that the exponent is sensitive to the system 
set-up, namely to its past history. This behaviour has already been observed in 
others numerical experiments at the infinite bifurcation point [F. Baldovin and 
A. Roblcdo, 2002b; U. Tirnakli et al., 1999]. 

A main difference is observed in the case of the tangent bifurcation. We 
repeated the same experiment for a = a tg with all the N initial conditions 
arranged in a single cell and let them expand entropically until the process of 
convergence to the attractor begins. After a transient time we found exactly the 
same power-law index as in the case of initial conditions uniformly distributed 
over the phase space (Fig. 5b). In this case the system seems to be not sensitive 
to the system set-up or to its past history. 

3 Conclusions 

In this paper we numerically studied how a statistical ensemble of initial con- 
ditions converges to the critical attractor using the logistic map. We compared 
the behaviour in two selected critical points: and a tg - The main results can 
be summarised in the following points: 

- for a = aoo a power-law regime with superimposed log-log periodic oscillations 
sets up after different crossover times, depending on the sampling ratio r and 
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there is no strict numerical evidence for the power-law exponent to be a constant 
independent on the sampling ratio, at least for r < 100. 

- for a = a tg a power-law convergence to the critical attractor, with different 
exponent, appears after a transient time even if the attractor has not fractal 
structure. The crossover time depends on the sampling ratio r in the same fash- 
ion as in the and is partition independent; 

- the choice of an ensemble of initial conditions concentrated on a single cell 
presents the same power-law exponent in the a tg case, whilst a different regime 
of convergence to the critical attractor appears in the case a = a^. 
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Figure 1 

We report for a — , in a log- log scale, the time evolution of the ratio between 
the occupied boxes and the available boxes. In frame (a) we show for different 
partitions W box of the interval [-1,1] the superimposed curves corresponding to 
the same r. We used the values W box = [20, 64, 128] x 10 3 for r = 0.1, W box = [5, 
10, 20, 64, 128] x 10 3 for r = 10 and W box = [10, 20, 128] x 10 3 for r = 1000. It is 
evident the occurrence of a power-law regime. Superimposed log-log oscillations 
appear for higher sampling ratios. Crossover times are also indicated. To the 
same sampling ratio r corresponds the same power-law exponent for the various 
grids until saturation is reached at different times. In frame (b) we show in 
sequence the values r = 0.1, 1, 10, 100, 1000 and 10000 for the two grids W box 
= 10000 (red) and W box = 128000 (black). The indicated crossover times show 
a net dependence on the sampling ratio r. 
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Figure 2 

We show the ratio W(t)/t _a71 versus log(t) for the infinite bifurcation point. 
In frame (a) we used W box = [10 (black), 20 (blue), 128 (rcd)]x 10 3 for different 
sampling ratios r. Frame (b) shows different sampling ratios calculated using 
Wbox = 20000 (red) and Wbox = 128000 (black). Higher sampling ratios exhibit a 
decreasing mean value of W(t)/t -0 - 71 , while the lower ones exhibit an increasing 
mean value, suggesting a slight dependence of the power-law exponent on r. 
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Figure 3 

We report for a — a tg , in a, log- log scale, the time evolution of the ratio of occu- 
pied boxes to the available boxes. In frame (a) we show for different partitions 
of the interval [-1,1] the superimposed curves corresponding to the same r. We 
used the values W box = [20, 64, 128] x 10 3 for r = 0.1, W box = [5, 10, 20, 64, 
128] x 10 3 for r= 10 and W box = [10, 20, 128] x 10 3 for r = 1000, as in the aoo 
case. It is evident the occurrence of a power-law regime well defined for higher 
sampling ratios with not superimposed log-log oscillations. Crossover times are 
indicated. Again, to the same sampling ratios r corresponds the same power-law 
exponent for the various grids. In frame (b) we show in sequence the values r = 
0.1, 1, 10, 100, 1000 and 10000 for the two grids W box = 10000 (red) and W box 
= 128000 (black). The indicated crossover times show again a net dependence 
on the sampling ratio r. The power-law regime is well pronounced at higher r 
values and extends in time over different orders of magnitude. Saturation effects 
appear for lower grids. 
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Figure 4 

Frame (a) evidences the power-law regimes for the tangent bifurcation point, 
with W box = [10 (blue), 128 (black)] x 10 3 and r = 1000, and with W box = 
10000 and r = 10000 (blue). Red lines are fits to the curves, with exponents /U 
= -1.55 for r = 1000 and fi = -1.45 for r = 10000. Different exponents describe 
the power-law regime for different sampling ratios. The power-law persists long 
in time over different orders of magnitude. Frame (b) shows the evolution of 
W(t)/t~ 155 versus log(t) for various r and grids W box = 128000 (black) and 
Wbox = 10000 (other colours). The exponent fi = -1.55 is appropriate only for 
the case r = 1000 where W(t)/t -155 shows a constant mean value. 
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Figure 5 

We compare the different behaviours for a — (frame (a)) and a = cit g (frame 
(b)), when all the N initial conditions are randomly selected inside a single bin 
of the phase space partition. We used in both cases Wbox = 64000. Frame (a) 
illustrates the two different power-law coefficients obtained for a = when the 
initial conditions are uniformly spread all over the phase space or are constrained 
inside a single bin. In the frame we show four different selections of the initial 
bin in different colours. The four curves overlap in the regime of convergence 
to the attractor, oscillating at high frequency among two states every temporal 
step and giving rise to the broad coloured curves. In frame (b) the two different 
selections of initial conditions exhibit the same power-law exponent. 
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Figure 1: We report for a — ctoo, in a log- log scale, the time evolution of the 
ratio between the occupied boxes and the available boxes. In frame (a) we 
show for different partitions Wf, ox of the interval [-1,1] the superimposed curves 
corresponding to the same r. We useS^the values Wbox = [20, 64, 128] x 10 3 
for r = 0.1, W box = [5, 10, 20, 64, 128] x 10 3 for r = 10 and W box = [10, 20, 
128] x 10 3 for r = 1000. It is evident the occurrence of a power-law regime. 
Superimposed log- log oscillations appear for higher sampling ratios. Crossover 




Figure 2: We show the ratio W(t)/t~ - 71 versus log(t) for the infinite bifurcation 
point. In frame (a) we used W box — [10 (black), 20 (blue), 128 (red)]x 10 3 for 
different sampling ratios r. Frame (b) shows different sampling ratios calculated 
using Wbox = 20000 (red) and W box 28000 (black). Higher sampling ratios 
exhibit a decreasing mean value of W(t)/t -0 ' 71 , while the lower ones exhibit an 
increasing mean value, suggesting a slight dependence of the power-law exponent 
on r. 




Figure 3: We report for a = at g , in a log- log scale, the time evolution of 
the ratio of occupied boxes to the available boxes. In frame (a) we show for 
different partitions of the interval [-1,1] the superimposed curves corresponding 
to the same r. We used the values Max = [20, 64, 128] x 10 3 for r = 0.1, 
W box = [5, 10, 20, 64, 128] x 10 3 for r = 10 and W box = [10, 20, 128] x 10 3 
for r = 1000, as in the Ooo case. It is evident the occurrence of a power-law 
regime well defined for higher sampling ratios with not superimposed log-log 



(a) (b) 




Figure 4: Frame (a) evidences the power-law regimes for the tangent bifurcation 
point, with W box = [10 (blue), 128 (black)] x 10 3 and r = 1000, and with W box 
= 10000 and r = 10000 (blue). Red lines are fits to the curves, with exponents [i 
= -1.55 for r — 1000 and fj, = -1.45 for^= 10000. Different exponents describe 
the power-law regime for different sampling ratios. The power-law persists long 
in time over different orders of magnitude. Frame (b) shows the evolution of 
W(t)/t~ 155 versus log(t) for various r and grids Wh nr — 128000 (black) and 




Figure 5: We compare the different behaviours for a = aoo (frame (a)) and 
a = a tg (frame (b)), when all the N initial conditions are randomly selected 
inside a single bin of the phase space partition. We used in both cases Wb ox 
= 64000. Frame (a) illustrates the two^fefferent power-law coefficients obtained 
for a — when the initial conditions are uniformly spread all over the phase 
space or are constrained inside a single bin. In the frame we show four different 
selections of the initial bin in different colours. The four curves overlap in the 



